Using the Standard Solar Model to Constrain Composition and S-Factors 
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While standard solar model (SSM) predictions depend on approximately 20 input parameters, 
SSM neutrino flux predictions are strongly correlated with a single model output parameter, the 
core temperature T c . Consequently, one can extract physics from solar neutrino flux measurements 
while minimizing the consequences of SSM uncertainties, by studying flux ratios with appropriate 
power-law weightings tuned to cancel this T c dependence. We re-examine an idea for constraining 
the primordial C+N content of the solar core from a ratio of CN-cycle 15 to pp-chain 8 B neutrino 
fluxes, showing that nonnuclear SSM uncertainties in the ratio are small and effectively governed 
by a single parameter, the diffusion coefficient. We point out that measurements of both CN-I cycle 
neutrino branches - ls O and 13 N /3-decay - could in principle lead to separate determinations of the 
core C and N abundances, due to out-of-equilibrium CN-cycle burning in the cooler outer layers of 
the solar core. Finally, we show that the strategy of constructing "minimum uncertainty" neutrino 
flux ratios can also test other properties of the SSM. In particular, we demonstrate that a weighted 
ratio of 7 Be and 8 B fluxes constrains a product of S-factors to the same precision currently possible 
with laboratory data. 

PACS numbers: 26.65.+t, 26.20.Cd,96.60.Fs, 14.60.Lm 



Important developments in solar neutrino physics have 
occurred over the past one to two years that impact the 
field's two major goals, probing the core of the Sun and 
constraining new weak-interaction phenomena. Supcr- 
Kamiokandc IV has reported a 8 B v flux measurement, 
(2.34 ± 0.03 ± 0.04) x 10 6 /cm 2 s [lj, that continues the 
progress toward high precision: combining in quadrature 
the statistical and systematic errors, one finds that the 
Super-Kamiokande uncertainty is now about 2%. The 
SNO combined analysis of solar neutrino data from all 
phases, has significantly narrowed the allowed range for 
the mixing angle #12 0], and new reactor and accelerator 
neutrino results have fixed the contributions of the sub- 
dominant mixing angle #13 Q . A new round of standard 
solar model (SSM) calculations has been completed to ex- 
plore competing compositions that optimize SSM agree- 
ment cither with the solar interior properties (as deter- 
mined from helioseismic mappings of the sound speed) or 
with solar surface properties (the interpretation of pho- 
toabsorption lines using our most sophisticated model 
of the Sun's atmosphere) The nuclear physics of 

the SSM has also been updated, with completion of the 
nuclear astrophysics community's second decadal evalu- 
ation of SSM S-factors and weak interaction rates [Ej]. 
Finally, Borexino has produced a precise (4.5%) mea- 
surement of the 862 keV neutrinos from 7 Be decay and 
first results on the pep flux, following its successful cali- 
bration campaign [6j. The 7 Be result has sharpened the 
"luminosity constraint" on the pp/pep neutrinos, which 
currently provides our most precise constraint on these 
fluxes if one assumes a steady-state Sun. 



The SSM, despite its relatively simple underlying 
physics, depends on ~ 20 input parameters, including 
the solar age and luminosity, the opacity, the rate of dif- 
fusion, the zero-age abundances of key elements (He, C, 
N, O, Ne, Mg, Si, S, Ar, Fe), and the S-factors for the 
pp chain (responsible for 99% of solar energy generation) 
and CN cycle. While SSM neutrino flux predictions gen- 
erally change with variation in any of these input pa- 
rameters, it has been recognized for some time that flux 
predictions are strongly correlated with a single output 
parameter, the core temperature T c Q. That is, a multi- 
dimensional set of variations in SSM input parameters 
{A/3j} from the SSM best values {/3f SM } often collapses 
to a one-dimensional dependence on T c , where T c is an 
implicit function of the variations {A/3,}: the effect of 
any variation {A(3j} on <f>i can be estimated simply from 
its effects on T c . The dominance of T c as the control- 
ling parameter for neutrino fluxes reflects the sensitiv- 
ity of Maxwellian-averaged rates (va(E)} to temperature. 
Consequently, if ^ oc T c 11 and 4>i 2 oc T c ' 2 , one can form 

weighted ratios 4>ii X ' 2 that are nearly independent 
of T c and thus nearly independent of variations over the 
multi-dimensional parameter space {A/3.,}. 

The situation becomes more interesting when two 
fluxes, say 4>ii and 4>i 2 i m addition to their common de- 
pendence on most underlying parameters {/3j}, have a 
very different dependence on some specific parameter, 
say j3\. In that case, from a weighted ratio of observed 
fluxes, one might be able to learn something about /3i, 
while the dependence on other input parameters largely 
cancels out. An example worked out previously Q and 
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updated below, is the additional linear dependence of CN 
neutrino fluxes on the primordial core number densities 
of C and N. To the extent that SSM uncertainties exceed 
the uncertainty of neutrino flux measurements, new in- 
formation can be extracted from neutrino measurements 
- in this case, the abundances of C and N in the Sun's 
primordial core. Another example we will discuss is the 
possibility that solar neutrino fluxes can be used to cross- 
check laboratory measurements of S-factors. 

The correlations between fa and T c are strong but 
not exact: for example, as the neutrino-producing core 
is extended (with the extent depending on the neutrino 
source), fluxes must depend on an integral over a core 
temperature profile, which cannot be exactly propor- 
tional to T c for all variations. Modern SSM calcula- 
tions allow one to address such issues, and to include 
their effects in the analysis. Monte Carlo studies can be 
done over wide classes of parameter variations {A/3j }, 
determining not only the best power-law descriptions of 
the fluxes, but to also the extent of reasonable variations 
around the power-law estimate. 

The sensitivity to parameter variations can be be ex- 
pressed in terms of the logarithmic partial derivatives 
a(i,j) evaluated for each neutrino flux fa and each SSM 
input parameter /3j, 



a(i,j) 



8ln[fa/$ SM ] 
dln[/y/3f M ] 



(1) 



where (/>f SM and /3| SM denote the SSM best values. This 
information, in combination with the assigned uncertain- 
ties in the /3j, then provides an estimate of the uncer- 
tainty in the SSM prediction of fa. Here we employ the 
logarithmic partial derivatives of which were evalu- 
ated for two different metallicities, corresponding to the 
higher Z composition of @, denoted GS98, and the lower 
Z composition of 0, denoted AGSS09. The older GS98 
abundances were obtained from a simple analysis of the 
solar atmosphere and yield excellent agreement with in- 
terior helioseismology. The newer AGSS09 abundances, 
obtained from a more sophisticated 3D model of the solar 
atmosphere that significantly improves the agreement be- 
tween measured and observed lines, are ~ 30% lower, and 
produce SSM sound speed profiles in significant disagree- 
ment with helioseismology. The SSMs evolved from these 
compositions are denoted SFII-GS98 and SFII-AGSS09 
in this paper, where SFII (Solar Fusion II) indicates the 
use of the latest nuclear S-factors @ . 

The logarithmic partial derivatives for the SFII-GS98 
SSM are given in Tables HI and HI1 divided as in H into two 
sets, corresponding to nuclear and "environmental" /3j-s. 
The nuclear parameters are the S-factors for the pp chain 
and CN cycle: Sn (p+p /3 decay), S 33 ( 3 He( 3 He,pp) 4 He), 
S 34 ( 3 He( 4 He, 7 ) 7 Be), S17 ( 7 Be( P , 7 ) 8 B), S e7 ( 7 Be elec- 
tron capture), and S114 ( 14 N(p,7) 15 0). The "environ- 
mental" parameters are those that directly influence the 



local temperature in the Sun, e.g., through their effects 
on evolution, the opacity, or SSM boundary conditions. 
They include the luminosity L Q , the Sun's age, the dif- 
fusion parameter, and the opacity. They also include the 
mass fractions of the principal solar metals, C, N, O, Ne, 
Mg, Si, S, Ar, and Fe, which have a significant influence 
on the opacity due to the strength of free -f-> bound tran- 
sitions. 

The partial derivatives allow one to define the power- 
law dependencies of neutrino fluxes, relative to the SSM 
best- value prediction </>? 



J.SSM 



N 

n 



where 



_ ^3 



Pi 



(2) 



and where the product extends over 19 SSM input pa- 
rameters of Tables [J and HH These derivatives determine 
how SSM flux predictions will vary, relative to 



iSSM 



, as 



the f3j are varied from their SSM best values. 

Our first example of the use of the logarithmic partial 
derivatives follows Q, though the analysis here differs 
in certain respects. The Sun produces about 1% of its 
energy through the CN-I cycle, which produces vs from 
the reactions 13 N(/3+) 13 C and 15 0(/3+) 15 N, with respec- 
tive /3-decay endpoints of 1.20 MeV and 1.73 MeV. Their 
fluxes in the SFII-GS98 SSM are 



(/>( 13 N) = 2.96(1 ± 0.14) • 10 8 /cm 2 s 
fa 15 0) = 2.23(1 ± 0.15) • 10 8 /cm 2 s. 



(3) 



The primordial C and N in the solar core are the cata- 
lysts for the conversion of four protons to 4 Hc via the CN 
cycle: the CN cycle alters the ratio of C to N as its burns 
into equilibrium, but does not change the total number 
density of C+N. The additional linear dependence of the 
CN-cycle on metallicity, due to this dependence on pri- 
mordial C and N, can be isolated by forming a ratio of 
fluxes that is effectively independent of T c , under varia- 
tions in all other SSM parameters. The appropriate ratio 
can be identified either by SSM Monte Carlo studies, at 
considerable cost numerically, or estimated from the log- 
arithmic partial derivatives. It was shown in Q that the 
two approaches yield essentially the same answer. 

The neutrino flux ratio identified in this way has the 
requisite metal sensitivity to distinguish GS98 abun- 
dances from those of AGSS09, resolving the solar abun- 
dance problem. More fundamentally, it will allow us to 
make a important test of a key SSM assumption, that 
the primordial Sun was homogeneous when nuclear burn- 
ing began - an assumption not obviously correct, given 
what we have learned in the past decade about large- 
scale metal segregation in the protoplanetary disk Q , but 
nevertheless critical to the SSM, which uses solar surface 
abundances to fix core abundances in the primordial Sun. 
Several groups have discussed relaxation of this assump- 
tion as a possibility for reconciling helioseismic data with 
AGSS09 abundances [l[ll|-[l3|. 
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TABLE I: Partial derivatives a(i, j) of neutrino fluxes with respect to solar environmental parameters and S-factors. Table 
entries are the logarithmic partial derivatives a(i,j) of the solar neutrino fluxes <f>i with respect to the indicated solar model 
parameter ft, taken about the SFII-GS98 SSM best values 0]. Several flux ratios that reduce the solar environmental factors 
are shown. 





Environmental ft 


Nuclear ft 


Source 


L 


Opacity 


Age 


Diffusion 


Su 


S33 


S34 




Se.7 


0114 


0(PP) 


0.766 


-0.112 


-0.100 


-0.013 


0.105 


0.034 


-0.067 


0.000 


0.000 


-0.007 


0(pep) 


0.989 


-0.318 


-0.024 


-0.019 


-0.217 


0.049 


-0.097 


0.000 


0.000 


-0.010 


0( 7 Be) 


3.434 


1.210 


0.760 


0.126 


-1.024 


-0.428 


0.853 


0.000 


0.000 


-0.001 


0( 8 B) 


6.914 


2.611 


1.345 


0.267 


-2.651 


-0.405 


0.806 


1.000 


-1.000 


0.007 


0( 13 N) 


4.535 


1.487 


0.932 


0.337 


-2.166 


0.031 


-0.062 


0.000 


0.000 


0.747 


0( 15 O) 


5.942 


2.034 


1.364 


0.382 


-2.912 


0.024 


-0.052 


0.000 


0.000 


1.000 


0( 7 Be)M 8 B) a465 


0.219 


0.002 


0.135 


-0.004 


0.209 


-0.240 


0.478 


-0.465 


0.465 


-0.004 


t /»( 13 N)/0( 8 B) - 576 


0.553 


-0.017 


0.157 


0.183 


-0.639 


0.264 


-0.526 


-0.576 


0.576 


0.743 


<?!»( 15 0)/^( 8 B) - 785 


0.515 


-0.016 


0.308 


0.172 


-0.831 


0.342 


-0.685 


-0.785 


0.785 


0.995 


0( 13 N)M 15 O) ' 776 


-0.075 


-0.091 


-0.126 


0.041 


0.093 


0.012 


-0.022 


0.000 


0.000 


-0.029 



TABLE II: As in Tabled] but for the partial derivatives a(i,j) with respect to the fractional abundances of the primordial 
heavy elements. . 





C, N ft 


Environment Abundance ft 


Source 


C 


N 


O 


Ne 


Mg 


Si 


s 


Ar 


Fe 


<Kpp) 


-0.008 


-0.002 


-0.006 


-0.005 


-0.004 


-0.010 


-0.007 


-0.002 


-0.021 


</>(pep) 


-0.016 


-0.003 


-0.012 


-0.006 


-0.003 


-0.013 


-0.015 


-0.005 


-0.062 


<M 7 Be) 


0.002 


0.001 


0.062 


0.055 


0.050 


0.104 


0.076 


0.019 


0.207 


<M 8 B) 


0.027 


0.007 


0.139 


0.109 


0.092 


0.192 


0.140 


0.035 


0.502 


<M 13 N) 


0.856 


0.165 


0.082 


0.058 


0.049 


0.111 


0.081 


0.021 


0.294 


<« 15 o) 


0.815 


0.217 


0.112 


0.081 


0.069 


0.150 


0.109 


0.028 


0.397 


<M 7 Be)/<M 8 B) - 465 


-0.011 


-0.002 


-0.003 


0.004 


0.007 


0.015 


0.011 


0.003 


-0.026 


4>( 13 N)/«X 8 B) - 582 


0.840 


0.161 


0.002 


-0.005 


-0.004 


0.000 


0.000 


0.001 


0.005 


<« 15 0)M 8 B) - 785 


0.794 


0.212 


0.003 


-0.005 


-0.003 


-0.001 


-0.001 


0.001 


0.003 


^( 13 N)/0( 15 O) - 776 


0.224 


-0.003 


-0.005 


-0.005 


-0.005 


-0.005 


-0.004 


-0.001 


-0.014 



The CN-cycle neutrino fluxes can be used as a direct 
probe of the core C and N abundances only to the extent 
that other SSM uncertainties can be controlled. Uncer- 
tainties in S-factors can in principle be improved through 
better laboratory measurements. In contrast, there may 
be no effective strategy to reduce the "environmental" 
uncertainties. These uncertainties often primarily affect 
the core temperature. For example, when metal abun- 
dances are varied, the SSM core temperature responds 
to the resulting changes in opacity and mean molecu- 
lar weight: high mctallicity cores are hotter. Neutrino 
fluxes also respond, reflecting their underlying power-law 
dependences on temperature. 



The dependence of the fluxes on environmental and 
other parameters can be determined from the logarithmic 



derivatives of Tables Q] and ITTl 
c/>( 13 N) 



13KHSSM 



c/>( 13 N) 



|^4. 535q1.487^0. 932 £)0.337j 



v fq-2.166 q0.031 q-0.062 qO.O cO.O cO.747] 
x L°H 33 °34 °17 °e7 a 114 J 

856 0. 165 0. 082 0.058 0. 049 0.111 J.081 0. 021 n. 2941 



'N 



c 6 x Ne x Mg X S 



b Fc J 



(4) 

where each parameter on the right-hand side represents a 
Pj/f3j M . The luminosity, opacity, solar age, and the dif- 
fusion parameters are denoted by L@, O, A, and D, while 
S and x denote S-factor or abundance ratios. Similarly, 

0( 15 O) 



isn'iSSM 



<K 15 o) 



|^5. 942^2. 034^1. 364 j-)0. 382j 



X 



[q-2.912 q0.024 q-0.052 qO.O qO.O ql.001 

L°ll D 33 D 34 °17 D e7 °114 \ 



[^815^0.217^0.112^0.081^0.069^.0.150^0.109^,0.028^.0.3971 



N 



■r 



o 



X 



No 



X 



Mg ^Si x Ar -^Fe J • 

(5) 

The 15 vs are of more interest experimentally, because 
their higher energy provides a window for observation in 



x 
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a scintillation detector, as discussed in [8(. From these 
expressions and from the "reasonable ranges" for input 
SSM parameters given in Table IIII1 one can then iden- 
tify the principal sources of SSM uncertainty in neutrino 
flux predictions. The ranges assigned in Table Hill to the 
metal abundances are of particular concern: the large 
differences between GS98 and AGSS09 reflect the ten- 
sion between helioseismology and 3D modeling of the so- 
lar atmosphere. For each metal, we assign to its abun- 
dance an uncertainty formed by two contributions. On 
one hand, following [jtjj . a systematic component based 
on the differences between the GS98 and AGSS09 com- 
positions given by 



Aft 

ft 



Abundance? S98 - Abundancef GSS09 



Abundance 



GSfJ.s 



Abundancef GSS09 



(6) 



to which we add in quadrature the observational uncer- 
tainty taken as the uncertainty quoted in the latest solar 
abundance compilation [loj| . This is a conservative ap- 
proach for assigning abundance uncertainties but it is 
appropriate for the present work, as it leads to robust 
upper limits to the precision with which solar neutrino 
experiments can constrain solar interior properties. 

It is reasonable to treat the effects of abundances in 
Eqs. (0} and as an overall scaling of mctallicity, as the 
differences between the GS98 and AGSS09 abundances 
are effectively systematic in the net mctallicity With 
this assumption, the dominant SSM uncertainty in Eq. 
([5]) is the core abundance of C + N which, if changed sys- 
tematically over a range equivalent to the GS98-AGSS09 
difference, alters the 15 neutrino flux by 30.7%. This 
is the sensitivity we want to exploit, in using CN neu- 
trinos as a probe of core metallicity. The next largest 
certainty comes from the If environmental parameters, 
16.5%: thus the environmental uncertainties are the pri- 
mary factor inhibiting our use of neutrinos as a probe of 
composition. The uncertainty coming from the S-factors, 
7.7%, is entirely dominated by S114, which alone con- 
tributes 7.2%. 

Now the nuclear physics uncertainties can be reduced 
with effort: in Q possible steps to improve existing mea- 
surements of S114 are described. But we have less con- 
trol over the environmental parameters, so an alternative 
strategy is needed to address these uncertainties. We use 
the well-measured flux of 8 B as a solar thermometer, to 
remove as much of the environmental dependence as pos- 
sible. 

We form a weighted ratio of the 15 O v and 8 B v fluxes 
to eliminate the dependence on T c to the extent possi- 
ble, and thus to minimize the dependence on 10 of the 
If environmental parameters of Tables HI and HH we do 
not include the diffusion coefficient, as this parameter 
plays a special role in the relationship between contempo- 
rary flux measurements and the primordial abundances 
we seek to constrain. The CN neutrino fluxes are more 



sensitive to diffusion, as Table U shows. All neutrino 
fluxes respond similarly to changes in core temperature 
induced by gravitational settling. However, the 15 O flux 
has an additional dependence on changes in the 12 C and 
14 N core abundances, as the rate is proportional to those 
abundances, for constant temperature. Thus the analysis 
is done in a way that isolates this additional dependence. 

To exploit the well-measured flux of 8 B neutrinos as a 
thermometer in this way, one must determine the linear 
correlations between ln(0( 13 N)) and ln(</>( 8 B)) and be- 
tween ln(</>( 15 0)) and ln((/>( 8 B)). While this can be done 
by direct Monte Carlo SSM calculations (see discussion 
below) , it was shown in Q that such an exercise is largely 
equivalent to minimizing the dependence on net logarith- 
mic derivatives. The solution to the minimization is a 
power law of the N observables, J|^ 1 ( ^f g i M ) b ' , with ex- 
ponents b\ given by the eigenvector with minimum eigen- 



value of the nuisance parameters error matrix [17 1 
( \ rvfi iWZ i\ 



Mi 



The computation of the matrices AWsb^n and A4sb. 15 o 
is straightforward. The direction of the smallest eigen- 
value is, respectively, 



</>("N) 



0(13 N )SSM 

x [S n ° 



&SSM(8 B ) 



0.576 



(3. 84(1 0.161 n 0.183 

Xq 



[L e 553 O- 



0.017 /10.1571 



X [Lq^U'"-"^ A" 

-0.639 n0.264 c-0.526 o -0.576 q0.576 
O33 3, 



J 34 



q-0.576 qO.576 cO.7431 
°17 D e7 D 114 J 



X 



0.002 -0.005 -0.004 O.0J. 0^,0. 001^,0.005 



o 



X 



and 

0( 15 O) 
0( 15 O) SSM / 



No 



x 



Mo 



Xg^ X<^ 



X 



Fc 



(8) 



hSSMfS 



( 8 b; 

X 



_ J. 794 0.212 n 0. 172 



n 0.785 

~C ,J} N 
515£)-0.016^0.308 







rc-0.831 nO.342 n-0.685 c-0.785 n0.785 nO.9951 
L°ll °33 °34 °17 a e7 °114 J 



0.003_-0.005 -0.003 -0.001 -0.001 „0. 001 _0.003 



'O 



x 



Nc 



X 



X 



Si 



x Ax x Fc 



(9) 



The important dependences on opacity and metallic- 
ity (other than C and N) have been almost entirely re- 
moved. While some residual dependence on luminosity 
and age remains, these parameters have relatively small 
uncertainties. Once the 8 B neutrino thermometer has 
removed the environmental effects, we find that the 15 O 
neutrino flux varies linearly under scaling of the 12 C and 
14 N abundances (0.794 + 0.212 = 1.006 - 1). This de- 
pendence can be made more explicit in Eq. (j9|) by the 
replacement 



,,0.794 0.212 



AS sm + N* SM 



(10) 
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TABLE III: Estimated la uncertainties in solar (from Bahcall, Serenelli, and Basu [3 and Fiorentini and Ricci and 
nuclear physics (from Adelberger et al. [f|) uncertainties, and their influence on flux predictions, computed from the partial 
derivatives of Table U 



ft 


Value 




Aci>( S B) ro/ x 

0(«b) y /o > 


Ac4( 7 Be) 
<K*Be) \ /o > 


<4(«N) \ /o > 


0(1S O ) ( /o > 


Lq 


3.842 xlO 33 ergs/s 


0.4 


2.8 


1.4 


1.8 


2.4 


Opacity 


1.0 


2.5 


6.5 


3.0 


3.7 


5.1 


Age 


4.57 Gyr 


0.44 


0.59 


0.33 


0.41 


0.60 


Diffusion 


1.0 


15.0 


4.0 


1.9 


5.1 


5.7 


P+P 


(4.01+0.04) xl0~ 25 MeV b 


1.0 


2.6 


1.0 


2.2 


2.9 


3 He+ 3 He 


(5.21±0.27) MeV b 


5.2 


2.1 


2.2 


0.16 


0.12 


3 He+ 4 He 


(0.56±0.03) MeV b 


5.4 


4.3 


4.6 


0.33 


0.28 


p+ 7 Be 


(20.8±1.6) eV b 


7.7 


7.7 


0.0 


0.0 


0.0 


e+ 7 Be 




2.0 


2.0 


0.0 


0.0 


0.0 


P+ 14 N 


(1.66±0.12) keV b 


7.5 


0.05 


0.0 


5.6 


7.5 



TABLE IV: Estimated la historical ("conservative") uncertainties in AGSS98 abundances, as defined in Bahcall and Serenelli 
[l6T | . The corresponding uncertainties in the neutrino fluxes are computed from the partial derivatives of Table [IT] 



ft 




<4(8b) y /o > 


A</>( 7 Be) (w \ 
<K»Be) { /0 > 


Ac4( 13 N) , W s 
<4(13N) V°) 


A0( 15 O) (o/s 
0(150) { /0 > 


c 


24.6 


0.66 


0.05 


21.1 


20.1 


N 


24.6 


0.17 


0.02 


4.1 


5.3 


O 


35.0 


4.9 


2.2 


2.9 


3.9 


Ne 


45.3 


4.9 


2.5 


2.6 


3.7 


Mg 


11.8 


1.1 


0.59 


0.58 


0.81 


Si 


11.8 


2.3 


1.2 


1.3 


1.8 


S 


13.8 


1.9 


1.0 


1.1 


1.5 


Ar 


34.9 


1.2 


0.66 


0.73 


0.98 


Fe 


11.8 


5.9 


2.4 


3.5 


4.7 



where Nq and Nn are the number densities of C and N. 
That is, the 15 neutrino flux depends effectively only on 
the sum of the number densities. The exponents appear- 
ing in Eq. ([9]) depend weakly on the SSM about which 
the variations are made. We use the SFII-GS98 SSM 
where N% SM /N^ SM ~ 0.80/0.20; the same ratio occurs 
in solar models using the solar composition from [loj . 

For 15 O, the case of most interest experimentally, the 
observable on the left hand side responds linearly to any 
scaling of the N and C primordial abundances. Diffu- 
sion, using Table Unl creates a 2.6% uncertainty in relat- 
ing contemporary flux measurements to the primordial 
abundance. This 2.6% is virtually all that remains of 
the original 16.5% SSM environmental uncertainty of Eq. 
(HJ): the 8 B neutrino thermometer has reduced the un- 
certainties associated with the remaining 10 parameters 
to below 0.35%. The third term on the right, contribu- 
tions from the S-factors, is now the dominant theoretical 
uncertainty in the relationship between primordial C+N 
and neutrino flux measurements, contributing 10.6% to 
the error budget. 

The explicit treatment of diffusion, effectively grouping 
diffusion with the C+N abundance, differs from the orig- 



inal work of [8[ . This choice is made for simple physical 
reasons, that neutrino flux measurements respond to con- 
temporary core abundances, yet the parameters needed 
in the SSM, which describes the Sun's evolution from 
the onset of nuclear burning, are primordial. Thus the 
relationship we establish between primordial core C+N 
and contemporary CN neutrino fluxes has a dependence 
on diffusion that should be made explicit, as we have 
done here. Indeed, the effects of diffusion are not incon- 
sequential: the SFII-GS98 and SFII-AGSS09 SSM metal 
profiles of Fig. [T]show that diffusion over 4.6 Gyr of so- 
lar evolution leads to nontrivial structures. Fortunately 
for our present goals, hclioseismology is sensitive to He 
and metal diffusion: the 15% uncertainty on the diffusion 
coefficient (see Table III) is a credible limit on diffusion 
uncertainties because of helioseismic constraints. 

Once this dependence on diffusion is separated out, it 
becomes apparent that almost all of the residual SSM 
"environmental" dependence identified in Q - variations 
in 11 SSM parameters producing a net uncertainty of 
2.6% - is due to the diffusion coefficient. The correla- 
tions illustrated in Fig. 3 of Ref. Q are redone in Fig. [3J 
with the removal of this one parameter. The net uncer- 
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tainty due to uncorrclatcd variations in the remaining 10 
parameters is now reduced to 0.3%, as mentioned above. 

The Super-Kamiokande measurement of the 8 B flux 
has reached a precision of 2%. Borexino has set the 
strongest constraint on the CNO solar neutrino inter- 
action rate (<7.9 counts/(day x 100 ton) at 95 % C.L.) 
and their latest purification campaign has resulted in a 
much lower background level, what opens the possibil- 
ity of the first detection of CNO neutrinos O]. SNO+ 
has the potential to measure the 15 O flux to an accuracy 
of about 10% in three years of running, if the detector 
design goals are reached Thus the current theo- 

retical ~ 30% uncertainty in the core C+N abundance 
could be substantially reduced by a neutrino measure- 
ment. In fact, the limiting uncertainty appears to be the 
nuclear physics, specifically S17 (7.7%) and S114 (7.5%). 
Both of these reactions were recently evaluated by the 
nuclear astrophysics community [H]. The uncertainty in 
S17 is dominated by the theory used to fit and extrap- 
olate measurements - the experimental contribution to 
the S17 error is 3.4% @. As ab initio methods may soon 
be available for such systems [2^, the situation could 
improve substantially. In the case of Si 14 a program of 
needed work was outlined in [f| , including new measure- 
ments to constrain the transitions to the 6.79 and 6.17 
MeV states in 15 O. We conclude that it should be pos- 
sible to significantly reduce the overall uncertainty from 
the nuclear physics. 

The expressions above are valid for the neutrino fluxes 
at the source. We need to account for the effects of neu- 
trino flavor conversion, as this alters the ratio of detected 
15 to 8 B neutrinos in detectors based on v x -e scattering. 
8 B neutrino oscillation probabilities are smaller because 
their energies correspond to the matter dominated flavor 
conversion while the CN neutrinos are in the vacuum os- 
cillation regime with small matter effects. We use the 
most up-to-date neutrino oscillation analysis [2~i| to es- 
timate the uncertainty due to the neutrino parameters, 
with lowers the weak-interactions uncertainty in our anal- 
ysis to ± 3%. 

The analysis suggests that a neutrino determination of 
the C + N content of the core at a confidence level of 
~ 10% is quite feasible with the future measurements. 
This assumes a 7% 15 O neutrino measurement and mod- 
est progress in lowering nuclear physics uncertainties to 
the same level. This should be compared to the cur- 
rent metallicity controversy, ~ 30%. Such a measure- 
ment would also constitute the first direct experimen- 
tal test of an important SSM assumption, that the pri- 
mordial core and modern solar atmosphere metallicitics 
are the same, once corrections are made for the effects 
of diffusion. Other sources of uncertainty - the Super- 
Kamiokande measurement of the 8 B neutrino rate for 
elastic scattering (ES), the SNO combined analysis con- 
straining weak interaction parameters, and the influence 
of diffusion (the one effect intrinsic to the SSM that can 



not be adequately subtracted using the 8 B neutrino ther- 
mometer) - are all of minor importance, contributing to 
the error budget at ~ 3%. 

There is a second constraint on core composition that 
could be obtained with CN neutrinos and that is inher- 
ently interesting because the observable is exceptionally 
free of SSM uncertainties. This constraint, however, re- 
quires a measurement of the 13 N neutrinos. If the pep 
shoulder is seen and if the level of background in the ES 
measurements can be kept low (or reliably subtracted), 
the remaining counts in the ~ 1 MeV region could be as- 
sociated with the 13 N and 15 neutrinos. One important 
observation is that the relative contributions of the two 
neutrino sources to the total rate vary with the electron 
recoil energy, due to the different neutrino energy distri- 
butions. This could allow the experimenters to separate 
the two CN-neutrino flux components. In Table [V] we 
show the relative contribution to the scattering rate as a 
function of the measured recoil electron energy interval. 
We have used the energy resolution of the Borexino detec- 
tor and the 13 N and 15 energy distribution of neutrinos 
given by the SFII-GS98 model. In the interesting energy 
range between the 7 Be and pep shoulders, the relative 
contribution to the rate varies by one order of magni- 
tude. The higher energy bins in this range are strongly 
dominated by the 15 O neutrino flux contribution and the 
error estimate of the flux will be comparable to the ex- 
perimental error in this energy region. The lower energy 
bins have a significant contribution from the 13 N neutrino 
flux. Thus this flux component could also be determined 
from the data, though with a larger uncertainty due to 
the need to separate this component from the dominant 
15 O contribution. An important caveat is the assumption 
there are no unidentified background sources that might 
mimic the 13 N neutrino signal. 

As noted previously, the CN cycle has not reached equi- 
librium in the Sun apart from its central core. The life- 
time of 14 N, determined by the rate of 14 N(p,7), is less 
than the solar age only for T7 > 1.33 (Tj is the tem- 
perature in 10 7 K). But at this temperature the lifetime 
of 12 C is ~ 2 • 10 7 years. Thus somewhat outside the 
central core, say at T7 ~ 1.15, there will be very lit- 
tle 14 N burning, and also very little 12 C burning, as the 
primordial carbon would have been consumed long ago. 
Still further outward, where T ~ 10 7 K, the 12 C lifetime 
is comparable to the solar age. This is the region in the 
contemporary Sun where primordial 12 C is being burned. 
We conclude that CN neutrinos are coming from two dis- 
tinct regions. The CN cycle is in equilibrium deep in the 
core, producing approximately equal numbers of 15 and 
13 N neutrinos, while well away from this region, in the 
cooler outer core at T ~ 10 7 K, primordial 12 C is burning 
to 14 N, producing only low-energy 13 N neutrinos. 

That is, the unequal fluxes of 13 N and 15 neutrinos 
are a reflection of the burning of primordial 12 C in the 
outer core. We can test for this effect by comparing these 
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FIG. 2: Solar neutrino fluxes from models in which all the parameters (black) or 10 of the 11 environmental parameters (red) 
- the diffusion coefficient is held fixed - are varied. At this resolution red points are indistinguishable from a line. The two 
upper panels show the correlation between the 8 B flux and the two CN-cycle neutrino fluxes. The slopes of the correlations 
between fluxes when only 10 environmental parameters are varied are given in the plots. The residuals from the fits, 0.31% 
and 0.32%, are shown in the lower panels. The results can be compared to those of Fig. 3 in @|, where residuals of 2.8% and 
2.6% were obtained for the 13 N and 15 fluxes, respectively, when diffusion was included as an 11th environmental parameter. 



8 



TABLE V: Ratio of the scattering rate by ls O and 1 N neutrino with electrons in a Borexino-like detector. 



E (MeV) 


[0.70,0.75] 


[0.75,0.80] 


[0.80,0.85] 


[0.85,0.90] 


[0.90,0.95] 


[0.95,1.0] 


i?( 15 0)/7?( 13 N) 


2.1 


2.6 


3.6 


5.5 


10.2 


24.7 



fluxes, treating the 15 O neutrino flux as the thermometer. 
In the exercise to find the linear correlation between the 
logarithmic fluxes, we now include diffusion among the 
environmental parameters, as it should affect 13 N and 
15 O neutrino rates almost equally. We find 
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(11) 



The residual environmental uncertainty is only ~ 0.7%. 
For the nuclear part we have made explicit the very small 
dependence on the S-factor of the 12 C(p,7) 13 N reaction 
rate, S112. The total nuclear uncertainty in the above 
expression is only 0.3%. We note here the dependence 
of this ratio on the nitrogen abundance is partially ac- 
cidental. The exponent that minimizes the environmen- 
tal uncertainties, 0.776, is very close to the ratio of the 
partial derivatives of these fluxes with respect to the N 
abundance 0.165/0.217=0.760 (see Table II). However, 
we find that even for unrealistically large variations of 
more than a factor of 2 in the assumed composition of 
the Sun, this ratio varies little, between 0.71 and 0.81. 
Therefore, the cancellation of the nitrogen abundance in 
Eq. [IT] will always occur at a level better than 0.5%. A 
similar conclusion can be drawn with respect to the Si 14 
astrophysical factor, which is also accidentally cancelled 
for the same reason. 

Other constraints: While we have focused on metal- 
licity and the CN-cycle neutrinos, due to the troubling 
solar abundance problem, the use of SSM power-law tem- 
perature dependences to extract parameter constraints 
is a general strategy for exploiting the Sun as a labora 
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tory. For example, the primordial 4 He abundance 
was recently constrained using very similar arguments. 
Another example we discuss here is the possibility of us- 
ing the SSM to cross-check laboratory measurements of 
S-factors. S17 is an important example because of the 
relatively large uncertainty in this S-factor and because 
of its importance to the branching between the ppll and 
ppIII cycles. The analysis is quite simple, a compari- 
son the the 7 Be and 8 B fluxes, neither of which has any 
anomalous dependence on metal abundances or diffusion. 
Thus we can optimize over all 13 non-nuclear parameters, 
with the anticipation that the residuals will be small in 



each of these. Following the previous calculation, we find 

0.465 



r Be) 



^Bc) 



SSM , 



4 SSM(8 B ) 



209 c-0.240 
D 33 



|^0. 219Q-0. 004^0. 135£)0.002j 

c-0.465 cO.465 q- 
a 17 °e7 a l 



qO.479 
D 34 



X 

0.0041 



X 



Mg ^Si 



Fc 



Sn 



i 0.24 r 



5:s:s 



S34S e 7 



'17 



0.48 



R 



nonnuc 



lcar 



SSM 



(12) 



The error introduced by grouping the astrophysical 
factors in the last expression is only 0.1%. In this ex- 
pression the factor Fgg™ uclcar represents the contribu- 
tions from the 13 non-nuclear uncertainties in the SSM: 
using the exponents above and the j3j of Tables III and 
IV, one finds that this contribution deviates from unity 
by ~ ±0.5%, and therefore plays no significant role. Thus 
effectively we have a direct relationship between neutrino 
flux measurements and nuclear cross sections. The left- 
hand side of Eq. (11) is the product of two factors. The 
first, [S11/S33] ' 24 , is uncertain to 1.3%, using the eval- 
uations of Solar Fusion II, with the error dominated by 
that in S33. The second, [S34S e 7/Si7] ' 48 , is uncertain 
to 4.6%, treating all uncertainties as uncorrelated. Thus 
the left-hand side of Eq. (11) is 1 ± 0.048, when all un- 
certainties are combined in quadrature. 

The right-hand side can be evaluated from the results 
of global solar neutrino flux analyses that incorporate the 
neutrino oscillation results important to mixing angle de- 
terminations (as the fluxes are the unoscillated instanta- 
neous ones) Q . The analysis is done in terms of the nor- 
malizations provided by SFII-GS98 SSM best values, for 
consistency with the logarithmic derivatives we employ: 
0( 7 Be) = 5.00 x 10 9 /cm 2 s and </>( 8 B) = 5.58 x 10 6 /cm 2 s. 
The experimental fluxes are 4.82(1 ± 0.045) x 10 9 /cm 2 s 
and 5.00(1 ± 0.03) x 10 6 /cm 2 s. Consequently the left- 
hand side is 1.016(1 ± 0.047). If the SFII-AGSS09 SSM 
best values are used to normalize the left hand side, the 
result is virtually unchanged, 1.015(1 ± 0.047). 

This result is significant: the constraint imposed on 



the ratio of S-factors [Sii/SjbP** [S 3 4S e7 /Si 7 ] - 48 rela- 
tive to SFII best values, using all information available 
from laboratory astrophysics, has the same precision as 
the similar ratio we can deduce from neutrino flux mea- 
surements, if we employ the SSM to predict the depen- 
dence of the fluxes on input S-factors, and if we con- 
strain all non-nuclear parameters to vary only within the 
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ranges allowed by their currently assigned uncertainties. 
This result was achieved by identifying a specific ratio 
of 7 Be and 8 B fluxes that the SSM predicts will exhibit 
the minimum uncertainty to variations in the 13 non- 
nuclear parameters. The two independent constraints - 
the left- and right-hand sides of Eq. (11) - are in excel- 
lent agreement, a result that reflects the concordance be- 
tween neutrino flux observations and laboratory nuclear 
cross section measurements, in the context of the SSM. 
As the uncertainty of the neutrino-flux result for this S- 
factor ratio is dominated almost entirely by that for the 
7 Be neutrino flux, further improvements in the Borcxino 
result (or new results from a next-generation experiment 
such as SNO+) would make the neutrino-flux S-factor 
constraint the more precise one. 

While this test of concordance between neutrino flux 
measurements and laboratory measurements of S-factors 
is our main point, one can be more aggressive and ask 
whether new S-factor information can be derived directly 
from neutrino flux measurements. As the most uncertain 
of the S-factors is S17, what level of precision is needed 
in neutrino flux measurements to improve our knowledge 
of this cross section? Equation (11) can be rewritten 



Si, 



S34S e 7 

0( 8 B) 



S11 



533 



0.5 



nonnuclcar 
SSM 



2.08 



ASSM/8 



(*B), 
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n 2.08 
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(13) 



Note that the simple exponent 0.5 on the S-factor ratio 
[S11/S33] is not accidental, but reflects the fact that the 
3 He abundance has achieved equilibrium in the region of 
the core where 7 Be and 8 B neutrinos are being produced. 
The number densities for 3 He and protons are then re- 
lated by 



A3 



cquil 



PP 



2A 



33 



(14) 



where X pp and A33 



are the local rates proportional to 
the respective S-factors. Effectively Eq. (fT3"| states 
that laboratory uncertainties in S17 (currently 7.5%) can 
be traded off against those in S34 (5.4%) and </>( 7 Be), 
the most poorly known quantities on the right-hand 
side. (Remember that all S-factors are normalized to 
their SFII best values.) Adding errors in quadrature, 
we find that this alternative determination yields S17 = 
0.967(1±0.117). The result will not be competitive, given 
the current laboratory precision of 7.5%, unless the un- 
certainties on both S34 and </>( 7 Be) are reduced to ~3%. 

Summary: We have refined the previous arguments of 
Q to show that future 15 O neutrino flux measurements 
have the potential to constrain the primordial core metal- 
licity of C+N to an accuracy of ~ 10%. This would be 
a very significant result, given that differences in recent 
abundance determinations exceed 30%. The method ex- 
ploits the additional linear dependence on metallicity of 



CN cycle burning, and is limited primarily by expected 
uncertainties of future experiments like SNO+ and by 
current uncertainties in laboratory measurements of nu- 
clear cross sections. The nonnuclear uncertainties in the 
relationship we derived, previously determined in Monte 
Carlo studies to be less than ~ 3% are in fact negli- 
gible apart from one parameter, the diffusion coefficient. 
The dependence on diffusion is natural, reflecting the fact 
that contemporary neutrino flux measurements are being 
used to constrain primordial abundances, not present-day 
core abundances. Our primary test of diffusion and its 
uncertainties comes from hclioscismology, which provided 
the initial motivation for including He and heavy-element 
diffusion in solar models. We also point out the possibil- 
ity - speculative experimentally, but intriguing theoret- 
ically - that by also measuring the 13 N solar neutrinos, 
one could determine the separate core abundances of C 
and N. The present-day burning of primordial C in the 
cooler outer core of the Sun contributes to the 13 N solar 
neutrino flux. 



The idea behind the metallicity extraction is a general 
one: forming ratios of neutrino fluxes that minimize the 
sensitivity to core temperature and thus to solar model 
uncertainties. We developed a second example of such a 
minimum-uncertainty SSM ratio - a comparison of 7 Be 
and scaled 8 B neutrino fluxes - that isolates a specific 
ratio of S-factors. We demonstrated that the precision 
to which this ratio is known from direct laboratory mea- 
surements is in fact identical to the precision it can be 
determined from measured solar neutrino fluxes and the 
SSM, given existing uncertainties on nonnuclcar input 
parameters to that model. Thus neutrino flux measure- 
ments have now reached the precision where meaningful 
consistency tests with laboratory cross sections can be 
done. In the example we developed, the laboratory cross- 
section measurements and neutrino fluxes were found to 
be in excellent agreement. 
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